1
1階微分方程式の構造解析
MATH009Lesson 2
00:00
現実のシステム、たとえば増加する貸付残高、落下する物体、あるいは絶滅危惧種の個体数を想像してください。 構造解析 1階微分方程式(ODE)の構造解析は、これらのシステムの将来の状態を予測するための数学的な橋渡しです。独立変数 $t$、従属変数 $y$、およびその瞬間的な変化率の関係を形式化します。

1. 構造的分類

本質的に、1階微分方程式は導関数と変数の関係を表します:$$\frac{dy}{dt} = f(t, y) \quad (1)$$ またはその陰関数形 $F(t, y) = 0$ です。方程式はその「骨格」によって分類されます。

  • 線形構造: 例えば $\frac{dy}{dt} = -ay + b$ (2) のような方程式で、関数が $y$ に関して線形である場合。注意:ここでは「一般解」という用語は、線形方程式にのみ適用します。
  • 自律的構造: 変化率が状態 $y$ にのみ依存する場合、$dy/dt = f(y)$ となります。このような方程式はしばしば 閾値レベル(T)という臨界個体数レベルを示すことが多く、これより下回ると種が拡散できず絶滅します。
  • 完全な構造: 条件 $M_y(x, y) = N_x(x, y)$ を満たすことで確認されます。この条件が成り立たない場合(例3のように)、系を満たす $\psi(x, y)$ は存在しません。
ステップ1:モデルの構築
物理的情報、たとえば 例4|脱出速度 (質量 $m$ の物体が地球から投げ出される)などは、数学的な表現に翻訳しなければなりません。重力と初期速度 $v_0$ を考慮しなければなりません。
ステップ2:安定性と存在性
私たちは リプシッツ条件: $|f(t, y_1) - f(t, y_2)| \le K|y_1 - y_2|$ によって、解の存在と一意性を保証します。この条件がないと、問題の「構造解析」が破綻したり、多価になる可能性があります。

2. 解と可視化

ある区間内のすべての $t$ に対して方程式を満たす微分可能な関数 $y = \phi(t)$ は、解と呼ばれます。幾何学的には、これを 積分曲線としてプロットします。ベルヌーイ方程式では、 変数変換 $v = y^{1-n}$ を用いて構造を線形化します。

🎯 重要な観察:オイラー法
例1( 例1 (年利12%の貸付残高 $S(t)$)において、オイラー法による離散近似 $y_{n+1} = y_n + f_n \cdot (t_{n+1} - t_n)$ は、実際の連続値よりも大きくなりがちです。これは解のグラフが 下に凸であるため、接線近似がグラフの上に位置してしまうからです。
$\frac{dS}{dt} = rS - k \implies y_n = \rho^n y_0$